Relaxation of Optically Excited Carriers in Graphene 
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We explore the relaxation of photo-excited graphene by solving a transient Boltzmann transport 
equation with electron-phonon (e-ph) and electron-electron (e-e) scattering. Simulations show that 
when the excited carriers are relaxed by e-ph scattering only, a population inversion can be achieved 
at energies determined by the photon energy. However, e-e scattering quickly thermalizes the carrier 
energy distributions washing out the negative optical conductivity peaks. The relaxation rates and 
carrier multiplication effects are presented as a function of photon energy and dielectric constant. 
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Following the success in fabricating field-effect tran- 
sistors using a single carbon atom layer [H, graphene 
has been actively studied as a possible basis for elec- 
tronic and optoelectronic devices [2]. Graphene's linear 
dispersion and zero energy gap result in unique trans- 
port and optical properties [3|-|5|, and new device con- 
cepts have been proposed to make use of those special 
features [gI, 0]- Recently, a new type of terahertz laser 
based on graphene has been proposed d, IH . Analytical 
model calculations considering relaxation involving op- 
tical phonon emission have shown that with sufficiently 
strong optical pumping, population inversion and neg- 
ative ac conductivity may be achieved in the terahertz 
frequency range This results from the fast optical 
phonon relaxation followed by an acoustic phonon re- 
laxation bottleneck for excited electrons and holes with 
energies E = Ep^Ti^o^ where Ep is the Fermi level, and 
Tiujo is the optical phonon energy. Some initial experi- 
mental results also have been reported suggesting stimu- 
lated terahertz emissions in graphene p^] . The threshold 
pumping power needed to obtain population inversion de- 
pends on the carrier recombination mechanism [s]. Pre- 
vious studies on quantum- well systems [11] have shown 
that electron-electron (e-e) scattering is responsible for 
the fast relaxation of photo-excited carriers, and recent 
experiments on graphene [HI, 13 1 also report ultrafast 
carrier kinetics. Therefore, it is essential to treat carrier 
dynamics rigorously to simulate the behavior of photo- 
excited carriers in graphene. 

Here we solve a time-dependent Boltzmann transport 
equation (BTE) for two-dimensional (2D) graphene with 
optical pumping. We consider carrier scattering and 
recombination-generation (R-G) mechanisms including 
acoustic and optical phonon scattering |14l4l6|. charged 
impurity scattering [17|] , photo-generation [4], sponta- 
neous optical transition [l8|, and e-e scattering that in- 
cludes impact ionization and Auger recombination (loj . 
Using the numerical solutions of the BTE, we study the 
influence of each scattering mechanism on the properties 
of optically pumped graphene, compare the results with 
those from previous approaches with phonon scattering 



only [s*], and address the challenges involved in realizing 
population inversion. We solve the BTE as 
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where f^^ is the distribution function at the nth time 
step for a state with 2D wavevector k = {kx^ ky) in the 
conduction {s = +1) or in the valence {s = —1) band, q 
is the unit charge, h is the reduced Planck constant, S'x 
is the electric field along the a;-direction, C/J^ ^ind G^^^ 
are the collision integral and the R-G rates, and At is 
the time interval, which is 0.05 ^ 10 fs in our simulation. 
We treat a bulk graphene under uniform illumination, so 
fs has no spatial dependence. 

The collision integral treats carrier scattering as 
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where ^i^k' is the transition rate from k to k', and 
ss' = ±1 means the intra- and inter-band transition re- 
spectively. For the two acoustic phonon modes, Fla and 
Fta, the sum of the two contributions can be treated as 
isotropic using an averaged sound velocity vs as 
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where Da is the acoustic phonon deformation potential, 
/3 is the 2D phonon wavevector with magnitude /3, pm 
is the mass density, A is the area of the sample, vp is 
the Fermi velocity, is the Bose-Einstein phonon oc- 
cupation number, S is the Dirac delta function, and ^ 
(or ±) signs represent the absorption and emission of 
phonons. For Da = 7.1 eV and vs = 17.3 km/s [16], we 
obtain a low field mobility of ~ 190,000 cm^V~^s~^ for 
carrier density of ~ 10^^ cm~^ as in Ref. Q. Acoustic 
phonon scattering is treated to be inelastic for a rigorous 
treatment of energy relaxation and broadening of dis- 
tribution functions. Electron-phonon (e-ph) coupling of 
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long- wavelengh optical phonon modes, Flo and Fto, is 
expressed as [IJ, ll5| 
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where Do = 3a/2Ji/2 ~ 11 eVA~^ is the optical phonon 
deformation potential with a coupling constant of Ji = 
5.3 eVA~^ [16], and TiuJo = 197 meV. For zone edge 
phonon modes, only the transverse mode Kto con- 
tributes to the carrier scattering, and the transition rate 
isgiven by Eq. (j2j) multiplied by (1 - ss' cos(^k — ^k0)/2 
[l^, with Do = 3Ji - 16 eVA-\ hujo = 157 meV, and 
6>k is the angle of k. 

While the charged impurity scattering [17] degrades 
the graphene mobility significantly due to momentum 
relaxation, we find that it has little effect on the opti- 
cal properties (the main focus of this work) because it 
does not relax energy. 

Coulomb scattering among electrons contributes to the 
R-G process through impact ionization and Auger re- 
combination [19], and it can also thermalize distribu- 
tion functions without changing the total carrier density. 
For Auger recombination and impact ionization, there 
are four processes involved, CCCV, VVCV, CVCC, and 
CVVV, where C and V represent the conduction and va- 
lence bands, and the second/third and first/fourth letters 
represent the initial and final states of the two electrons 
involved. For example, the transition rate from k to k' 
in the conduction band due to CVCC is expressed as 
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where \M\^ = |M^p + \M,\^ - \Md - M,\^ = 2\Md\^ + 
2 1 Me p - M* Me - M^M*, and M^ and Me are the direct 
and exchange matrix elements. To treat the K-K' val- 
ley degeneracy, we double the 2\Md\'^ and 2|Mep terms. 
Coulomb interactions are calculated using the random 
phase approximation (RPA) [l^. The resulting collision 
integral in Eq. ([1]) becomes four-dimensional, and it can 
be reduced to 2D using the momentum and energy con- 
servation conditions [19]. Processes conserving the num- 
ber of electron-hole pairs, i.e., CCCC, WW, CCW, 
and CVCV, are treated in a similar way. 

Gsk in Eq. ^ includes photo-generation and sponta- 
neous emission. The photo-generation rate is 
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where a is the fine structure constant, huj, 
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the incident photon energy, Popt is the pumping power, 
and ss' = —1. From Eq. (|3]), we can show that the 
graphene with f_^i=0 and /_i = 1 has the absorbance 
of na 2.3 % The spontaneous emission is treated 



using the approach in Ref. [18|, but its contribution is 
much smaller than those from other processes. 

Using the models above, we numerically solve the BTE 
on a 2D A:-grid. We first simulate optically pumped 
graphene without considering e-e scattering, where e-ph 
scattering is the main mechanism for energy relaxation as 
in the previous approaches [sl, ISi] . Then we introduce e-e 
scattering and see how the results change. We also ex- 
plore the effects of Eopt, Popt, and Ep on the properties of 
optically pumped graphene. All simulation results are for 
300 K, zero (ox^ and no charged impurity. Photo-current 
with finite (ox due to the transport of photo-excited car- 
riers will be discussed elsewhere. 

Figure [1] shows simulation results under steady-state 
low power optical pumping considering e-ph scattering 
only. In all following results, f±i means /±i,k averaged 
over k with |k| = /c and E = ±hvFk. In Fig. [T]^a) for 
Ep = eV, /+i under optical pumping dramatically devi- 
ates from equilibrium with peaks modulated by Eopt- For 
Eopt = 0.6 eV, for example, photo-excited electrons first 
generated at = 0.3 eV pile up at ~ 0.1 eV and 0.15 
eV due to the emission of Flq /to and Kto phonons. For 
Ep = eV, photo-excited holes pile up symmetrically in 
the valence band. Due to interband optical phonon scat- 
tering, the piled-up electrons and holes recombine with 
holes and electrons in the opposite bands, so the carri- 
ers with E < 0.1 eV are depleted in Fig. [U^a). Figure 
[Tfb) shows the optical absorbance due to interband tran- 
sition, 7ra(/_i — /+i), vs. frequency 2\E\/h, where h 
is the Planck constant. The absorbance becomes nega- 
tive, i.e., population inversion is achieved, at specific fre- 
quency ranges depending on Eopt- Simulation results for 
Ep = 0.05 eV and 0.1 eV are shown in Figs, mc)-(d). For 
finite Ep^s^ /+i and /_i are not symmetric, and the fre- 
quency and the height of the population inversion peaks 
depend on Ep. In general, as Ep increases, a smaller 
energy region becomes available for the minority carriers 
to experience an acoustic phonon bottleneck, so a larger 
pumping power is required to invert the population. 

Figure [2] shows simulation results for optically pumped 
graphene under steady-state with e-e scattering consid- 
ered. In Fig. [21(a) for Ep = eV and dielectric constant 
hi = 2.5, /+i with optical pumping does not deviate much 
from equilibrium, and there is little Eopt-dependence un- 
like in the case of Fig. [D^a). We find that under photo- 
excitation the distribution function remains equilibrium 
with an elevated electronic temperature determined 
from the Fermi-Dirac distribution. In Fig. Efa), = 381 
K, 402 K, and 419 K for Eopt = 0.4 eV, 0.6 eV, and 0.8 eV 
respectively. In Fig. EJb), although we increase Popt by 
five times and use a high = 10 to reduce the Coulomb 
interaction, the distribution functions are more broad- 
ened with higher Te's (497 K, 529 K, and 553 K) and 
still far from population inversion, and the absorbance 
never goes below zero in Fig. [2](c). Here we should note 
that in RPA [17^] , the Coulomb interaction does not scale 
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FIG. 1. Simulation results for optically pumped graphene 
under steady-state with e-ph scattering only, (a) /+i vs. E 
and (b) absorbance for Ef = eV, Eopt = 0.4/0.6/0.8 eV, and 



D^p, = 20/30/40 kWcm- 
for = 0.05/0.1 eV, Eopt 

Population inversion is achieved and modulated by E, 



(c) /+i vs. E and (d) absorbance 
0.8 eV, and Popt = 40 kWcm"^ 



opt ' 



as as expected from the unscreened modeL For Ef 
= eV, for example, the RPA coincides with the un- 
screened model with hz + q'^ /{ShvpSo) ^ + 3.5 where 
eo is the vacuum permittivity, so changing from 2.5 to 
10 results in the Coulomb interaction being reduced by 
about five times, not by 16 times. Simulation results for 
Ef = 0.05 eV and 0.1 eV are shown in Figs. [2ld)-(f). 
In Fig. [2jd), population inversion is never achieved for 
K = 2.5 and the same Popt^s as those in Fig. HJc), and 
the distribution functions are thermalized with Tg = 408 
K and 379 K for Ep = 0.05 eV and 0.1 eV respectively. 
With = 10 and larger pumping powers in Fig. [2fe)-(f), 
distribution functions are more broadened with elevated 
Tg's (536 K and 494 K) instead of showing multiple peaks 
as in Fig. Hl^c). 

Therefore, e-e scattering significantly broadens the dis- 
tribution functions washing out peaks due to cascade 
emission of optical phonons fu] . To observe this quanti- 
tatively, we calculate the carrier lifetime Tg^k = Xlk'(-'- ~ 
/s'kO*^k,k', which is the characteristic time it takes for a 
carrier in the state with k to scatter out to other states. 
Figure [3] shows results for vs. E for conduction band 
under equilibrium with Ef = 0.1 eV and 0.05 eV and 
= 10, where r is the average of Tk for |k| = /c. The con- 
tribution from e-e scattering scales linearly with E", and 
it is much larger than those from e-ph scattering and 
increases with decreasing Ef due to reduced screening. 

So far we have mainl y d iscussed steady- state results, 
but transient responses [20| are also important in under- 
standing carrier dynamics in graphene. In Figs. IUa)-(b), 
we compare the photocarrier density An = n — no vs. t 
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FIG. 2. Simulation results for optically pumped graphene 
under steady-state with e-e scattering considered, (a) /+i vs. 
E for = eV, At = 2.5, Eopt = 0.4/0.6/0.8 eV, and Popt 
= 20/30/40 kWcm"^ (b) /+i vs. ^and (c) absorbance with 
increased pumping powers and n — 10. (d) /+i vs. E'for Ef — 
0.05/0.1 eV, K = 2.5, Eopt = 0.8 eV, and Popt = 40 kWcm-^ 
(e) /+i vs. E and (f) absorbance with increased pumping 
powers and k — 10. Distribution functions are thermalized 
with elevated Te's, and population inversion is not achieved. 
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FIG. 3. Simulation results for vs. E (conduction band, k 
— 10) under equilibrium with (a) Ef = 0.1 eV and (b) Ef — 
0.05 eV. Coulomb scattering rates are much larger than e-ph 
scattering rates and increase with decreasing Ef- 



for Ef = eV with and without e-e scattering (k = 
2.5), where n is the total electron density, no is the equi- 
librium electron density, and the optical pumping starts 
at time t = ps. With e-e scattering, steady-state is 
quickly achieved, and An, which in turn determines the 
photo-current response [2^, is much smaller than with e- 
ph scattering only. We can fit the An vs. t curves using 
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FIG. 4. Simulation results for An vs. t with optical pumping 
(Eopt = 0.4/0.6/0.8 eV, Popt = 20/30/40 kWcm-^) for Ef 
= eV with (a) e-ph scattering only and (b) e-e scattering 
considered with At = 2.5. (c) Relaxation time vs. Eopt and 
(d) absorption efficiency vs. Eopt for various k^s for the same 
optical pumping conditions. 



In summary, we numerically solved a time-dependent 
BTE for bulk graphene considering photo-generation and 
relevant carrier scattering and R-G mechanisms to ex- 
plore properties of optically pumped graphene and pos- 
sibilities to achieve population inversion. Simulation re- 
sults showed that when carrier energies are relaxed by 
e-ph scattering only, it is possible to achieve population 
inversion at frequencies controlled by incident photon en- 
ergy and graphene Fermi level. A more realistic model 
considering Coulomb scattering, however, showed that 
distribution functions are significantly broadened, and it 
becomes hard to achieve population inversion even with 
much higher pumping power. The results stress the im- 
portance of e-e scattering in carrier dynamics in graphene 
and suggest that the Coulomb interaction among carri- 
ers must be significantly suppressed, such as by bandgap 
opening in bilayer graphene, to realize graphene-based 
terahertz lasers. However, e-e scattering leads to signifi- 
cant carrier multiplication. 
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sity for helpful discussions. Computational support was 
provided by the Network for Computational Nanotech- 
nology, supported by the National Science Foundation. 



An = Anrn{l — exp(— t/r^)) and extract the relaxation 
time Tn and the absorption efficiency r] = An/Ropt, where 
Ropt = Popt/ Eopt' For e-ph scattering only in Fig. H^a), 
Tn is very large (45 ~ 60 ps) and increases with Eopt 
while it is much shorter for e-e scattering (~ 3 ps) in 
Fig. m^b). As shown in Fig. HJ^c), increases with in- 
creasing and somewhat decreases with increasing Eopt- 
Figure m^d) shows the rj vs. Eopt results. For = oo (i.e. 
e-ph scattering only), 77 is close to na and independent of 
Eopt' For finite /^'s, however, r] is greater than 7ra, which 
indicates carrier multiplication. The rj increases with 
decreasing due to the increased Coulomb interaction, 
and it increases with increasing Eopt because high en- 
ergy photocarriers have higher ability to initiate impact 
ionization. A recent theoretical study has reported that 
with pulse excitations, impact ionization leads to dra- 
matically increased photocarrier densities [21]. However, 
for continuous excitation (the main focus of this work), 
e-e scattering leads to a lower steady-state An = AnTn 
because Auger recombination becomes stronger and re- 
combination due to optical phonon emission is boosted 
by the broadened distribution functions, which all result 
in a reduced r^. We find that Tg follows the similar time 
evolution as An, and the steady-state Tg under optical 
pumping depends on Eopt^ Popt^ and n as shown in Fig. 
[21 We also find that the decay time constant, after the 
optical pumping is switched off, is slightly different than 
the rise time constant due to the difference in the scat- 
tering rate determined by the distribution function in the 
photo-excited state. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 

[7] 

[8] 

[10: 
[11 

[12 
[13 
[14 
[15 
[16 

[17 

[18 

[19 
[20; 

[21] 



K. S. Novoselov et a/., Science 306, 666 (2004)" 
P. Avouris, Nano Lett. 1 0, 4285 (2010) 

Nature Nanotech. 3, 206 (2008). 
1308 (2008) 



J.-H. Ghen et al, 
R. R. Nair et al, 
K. F. Mak et al, 
M. I. Katsnelson 



Science 320, 

Phys. Rev. Lett. 101, 196405 (2008)] 
, K. S. Novoselov, and A. K. Geim, 
Nature Phys. 2, 620 (2006). 

V. V. Gheianov, V. Fal'ko, and B. L. Altshuler, 
[Science 315, 1252 (2007), 

V. Ryzhii, M. Ryzhii, and T. Otsuji, 

J. Appl. Ph ys. 101, 083114 (2007) 

V. Ryzhii et al,\J. Appl. Phys. 107, 054505 (2010) | 

H. Karasawa et al, J. Infrared Milli. Terahz. Waves , 1 

(2010). 

S. M. G oodnick and P. LugU, 

[Phys. Rev. 257 8 (1988) 

R. W. Newson et al, Opt. Exp ress 1 7, 2326 (2009) | 
H. Ghoi et a/.,Xpplied Physics Letters 94, 172102 (2009)] 
M. Lazzeri et al, Phys. Rev. Lett. 95, 236802 (2005) 
T. Ando, J. Phys. Soc. Jpn. 75, 124701 (2006). 

V. Perebeinos and P. Avouris, 

[Phys. Rev. B 81, 195442 (2010) [ and references therein. 

E. H. Hwang and S. Das Sarma, 
[Phys. Rev. B 77, 195412 (2008) 

F. T. Vasko and V. Ryzhii, 
|Phys. Rev. B~7, 195 43 3 (2008), 

F. Rana, |Phys. Rev. B 76, 155431 (2007)| 

T. Hertel and G. Moos, 

|Phys. Rev. Lett. 84, 5002 (2000) | 

T. Winzer, A. Knorr, and E. Malic, 

|Nano Lett. 1074839 (2010) 



